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RESUMEN 

El resumen sera traducido al espanol por los editores. The first detections of afterglows from short 
gamma-ray bursts (GRBs) have confirmed the previous suspicion that they are triggered by a different central 
engine than long bursts. In particular, the recent detections of short GRBs in galaxies without star formation 
lends support to the idea that an old stellar population is involved. Most prominent are mergers of either double 
neutron stars or of a neutron star with a stellar-mass black hole companion. Since the final identification of the 
central engine will only come from an integral view of several properties, we review the observable signatures 
that can be expected from both double neutron stars and neutron star black hole systems. We discuss the 
gravitational wave emission, the structure of the neutrino-cooled accretion disks, the resulting neutrino signal 
and possible mechanisms to launch a GRB. In addition, we address the speculative idea that in some cases a 
magnetar-like object may be the final outcome of a double neutron star merger. We also discuss possibilities 
to explain the late-time X-ray activity that has been observed in several bursts. 

ABSTRACT 

The first detections of afterglows from short gamma-ray bursts (GRBs) have confirmed the previous suspicion 
that they are triggered by a different central engine than long bursts. In particular, the recent detections 
of short GRBs in galaxies without star formation lends support to the idea that an old stellar population is 
involved. Most prominent are mergers of either double neutron stars or of a neutron star with a stellar-mass 
black hole companion. Since the final identification of the central engine will only come from an integral view 
of several properties, we review the observable signatures that can be expected from both double neutron 
stars and neutron star black hole systems. We discuss the gravitational wave emission, the structure of the 
neutrino-cooled accretion disks, the resulting neutrino signal and possible mechanisms to launch a GRB. In 
addition, we address the speculative idea that in some cases a magnetar-like object may be the final outcome 
of a double neutron star merger. We also discuss possibilities to explain the late-time X-ray activity that has 
been observed in several bursts. 
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Gamma-ray bursts (GRBs) are subdivided on 
grounds of both their durations and spectra into 
short-hard and long-soft bursts (Kouveliotou et al. 
1993). While a lot has been learned from the de- 
tection of long GRB afterglows since 1997 (Costa et 
al. 1997, van Paradijs et al. 1997, Frail et al. 1997, 
Bloom et al. 1999, Stanek et al. 2003, Hjorth et 
al. 2003), for example that long GRBs are related 
to death of massive stars in star-forming regions, the 
afterglows of short GRBs remained elusive until the 
early summer of 2005. The recent detections of af- 
terglows from short GRBs point to a central engine 



that is related to an old stellar population. In partic- 
ular, it has been argued (Bloom et al. 2006, Berger 
et al. 2005, Barthelmy et al. 2005, Villasenor et 
al. 2005) that the observations are consistent with 
being the result of compact binary mergers: they 
occur systematically at lower redshift^] than their 
long-duration cousins (e.g. Fox et al. 2005), both 
in galaxies with (Hjorth et al. 2005) and without 
star formation (Berger et al. 2005, Barthelmy et al. 
2005) and they are not accompanied by a detectable 
supernova explosion (Bloom et al. 2006, Fox et al. 
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2 It has recently been suggested (Berger et al. 2006b) that 
at least 1/4 of the short bursts could lie at redshifts z > 
0.7 and would therefore imply substantially larger isotropised 
energies than was inferred for the first set of detected short 
GRBs with afterglows. 
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Fig. 1. Merger of a double neutron star system with 1.4 
M© each (snapshots at 4.41, 5.54 and 15.0 ms). The low- 
est panel shows the gravitational wave amplitudes times 
the distance to the source, d (in units of 1.5 km). 

2005, Hjorth et al. 2005). 
There are, however, puzzling observations of late- 



time flaring activity in several short bursts Berger et 
al. 2006a). For example, GRB050724 showed long- 
lasting (~ 100s) X-ray flaring activity after a delay 
of ~ 30 s. These long-lasting flares motivated Mac- 
Fadyen et al. (2006) to suggest an accretion-induced 
collapse of a neutron star in a low-mass X-ray binary 
as an alternative central engine. This central engine 
would produce the flaring by interaction of the GRB 
outflow with the non-compact companion star. 
As the final identification of the central engine will 
probably only come from an integral view of various 
features from a large sample of events, we will give 
an overview over the properties and signatures that 
can be expected from compact binary mergers. They 
will include gravitational waves, the structure of the 
forming accretion disk, the neutrino emission and the 
burst itself. Moreover, we will discuss possibilities to 
produce the observed late-time flaring activity. 
The discussion is mostly based on own simulations, 
some of which have been reported previously (Ross- 
wog and Davies 2002, Rosswog and Ramirez-Ruiz 
2002 and 2003, Rosswog and Liebendorfer 2003, 
Rosswog et al. 2003, Rosswog et al. 2004, Ross- 
wog 2005) and several new simulations that are in- 
troduced here, see Table [U The physics included in 
our models and the numerical techniques that we use 
have been described in a series of papers (Rosswog 
and Davies 2002, Rosswog and Liebendorfer 2003, 
Rosswog et al. 2003, Rosswog et al. 2004). 

The new neutron star merger calculations have 
been produced with a new Smoothed Particle Hy- 
drodynamics (SPH) code that benefits from a slew 
of numerical improvements. It incorporates 

• an artificial viscosity oriented at Riemann 
solvers, see Chow and Monaghan (1997), 

• a consistent accounting of the effects of the so- 
called "grad-h" -terms, see Springel and Hern- 
quist (2002), Monaghan (2002) and Price 
(2004), 

• a consistent implementation of adaptive gravi- 
tational softening lengths, see Price and Mon- 
aghan (2006), 

• the option to evolve magnetic fields together 
with the fluid, either via so-called Euler poten- 
tials (Stern 1970) or via a SPH discretisation 
of the MHD-equations (Price and Monaghan 
2005). 

First results obtained with this code have been pub- 
lished (Price and Rosswog 2006), a detailed code 
documentation can be found in Rosswog and Price 
(2006). 
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TABLE 1 

NEW RUNS: DNS STANDS FOR DOUBLE NEUTRON STAR SYSTEMS, NSBH FOR NEUTRON STAR 
BLACK HOLE SYSTEMS. M 1 AND M 2 ARE THE COMPONENT MASSES, Q IS THE MASS RATIO, A 
IS THE INITIAL SEPARATION, N STANDS FOR NEWTONIAN, PW FOR THE Paczynski -WIITA 
PSEUDO-POTENTIAL, R ABS IS THE RADIUS OF THE ABSORBING BOUNDARY, IN UNITS OF 
GMbh/C 2 , THE NEXT COLUMN GIVES THE SPH PARTICLE NUMBER. 

run type Mi, M2, q ao [km] grav. R a b s SPH part. 



DA 


DNS 


1.4, 1.4 , 1 


48 


N 


DB 


DNS 


1.1, 1.6 , 0.6875 


48 


N 


NA 


NSBH 


1.4, 3.0 , 0.4667 


60 


N 


NB 


NSBH 


1.4, 4.0 , 0.35 


72 


N 


NC 


NSBH 


1.4, 7.0 , 0.2 


90 


N 


ND 


NSBH 


1.4, 10.0 , 0.14 


112.5 


PW 



2. GRAVITATIONAL WAVES AS PROBES OF 
THE COALESCENCE DYNAMICS 

Gravitational waves can serve as a direct probe 
of the merger dynamics and, if detected coincident 
with a short GRB, they would provide the ultimate 
proof of the compact binary nature of the central en- 
gine. The dynamics of the merger process is sensitive 
to the mass ratio of the involved components, there- 
fore, the signal from the merger of a double neutron 
star system (DNS) with a mass ratio close to unity 
can be very different from a neutron star black hole 
(NSBH), where the black hole can, in principle, be 
much more massive than the neutron star. 
In the last minutes before the coalescence the grav- 
itational wave signal will slowly sweep through the 
frequency range that is accessible to ground-based 
gravitational wave-detectors such as GEO600 (Grote 
et al. 2005), LIGO (Abramovici et al. 1992), 
TAMA (Ando et al. 2004) or VIRGO (Spallici et 
al. 2005, Freise 2005). The detection of a gravi- 
tational wave signal coincident with a short GRB 
would allow to determine the distance to the source, 
the luminosity and the beaming angle (Kobayashi 
and Meszaros 2003). 

The coalescence of a double neutron star system can 
be divided into three phases: the secular inspiral due 
to gravitational wave emission, the actual merger 
and, finally, the "ring-down" phase, in which the 
death-struggle of the freshly formed super-massive 
neutron star will take place. After this phase the 
remnant will have settled into its final state, proba- 
bly a black hole (an alternative is discussed in Sec. [6]). 
During the inspiral, the gravitational wave forms can 
be accurately described via post-Newtonian expan- 
sions for point masses. To date, post-Newtonian for- 
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malisms exist up to 3.5PN order, see Blanchet (2006) 
for a recent review, and lowest order spin-spin and 
spin-orbit couplings can be accounted for (e.g. Will 
2005). In the inspiral phase, both the frequency and 
the amplitude of the waves will increase, the system 
is said to "chirp". The orbit decays secularly until 
the last stable orbit is reached, where the binary en- 
ters a "plunging phase" . In this phase the stars fall 
nearly radially towards each other and merge within 
about one orbital period. This dynamical instability 
(Chandrasekhar 1975, Tassoul 1975) is the result a 
steepening of the gravitational potential due to both 
purely Newtonian tidal (Lai et al. 1993) and general 
relativist ic effects, for a recent review on relativis- 
tic binaries see Baumgarte and Shapiro (2003). For 
the last two phases, the merger and ringdown, three- 
dimensional hydrodynamic simulations are required 
to predict the gravitational wave signal. 
Double neutron star merger wave forms has been pre- 
dicted by several groups. The calculations started 
with Newtonian calculations (e.g. Oohara and Naka- 
mura 1989, Rasio and Shapiro 1994, Ruffert et al. 
1997, 2001). Post-Newtonian approaches (Ayal et 
al. 2001, Faber and Rasio 2000) turned out not to 
be particularly successful due to the importance of 
higher order PN-corrections. In the conformal flat- 
ness approach (Isenberg 1978, Wilson et al. 1996) 
it is assumed that the dynamical degrees of freedom 
of the gravitational fields, i.e. gravitational waves, 
can be neglected and that the spatial part of the 
metric is (up to a conformal factor) fiat and remains 
so during the further evolution. This approxima- 
tion has been used in several compact binary sim- 
ulations (e.g. Wilson et al. 1996, Oechslin et al. 
2002, 2004, Faber et al. 2004, Oechslin 2006) and 
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Fig. 2. DNS merger with 1.1 and 1.6 M (at 1.51, 6.17 
and 15.0 ms). The lowest panel shows the gravitational 
wave amplitudes times the distance to the source, d (in 
units of 1.5 km), as a function of time. 
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Fig. 3. Merger of 1.1 and 1.6 Mq system (t=1.63 ms): 
Impact of accretion stream onto surface of accreting neu- 
tron star. Color-coded is density in the orbital plane, the 
arrows indicate velocities. 

it could be shown (Cook et al. 1996) to be provide 
reasonable accuracy in several cases. Finally, fully 
relativistic simulations have been performed by Shi- 
bata and collaborators (Shibata 1999, Shibata and 
Uryu 2002, Shibata et al. 2005). For a recent review 
on numerical relativity and compact binaries we re- 
fer to Baumgarte and Shapiro (2003). 
The simulation of neutron star black hole binaries 
has been somewhat lagging behind the DNS status. 
Again, first simulations were Newtonian (e.g. Lee 
1999a,b, Janka 1999, Lee 2000, 2001, Rosswog et al. 
2004), slightly later followed simulations that used a 
Paczyhski -Wiita (1980) pseudo-potential (Lee and 
Kluzniak 1999a, Rosswog 2005, Setiavan et al. 2004, 
2006). Very recently, progress has been made with 
relativistic approaches. For example, Taniguchi et 
al. (2005) were able to construct general relativistic 
quasi-equilibrium sequences for black holes that are 
much more massive than the neutron star. Faber 
et al. (2006) used an exact metric together with a 
conformal flatness approach in which the black hole 
position was artifically kept fixed in space. Loffler et 
al. (2006) were able to treat binary components of 
comparable mass, but their approach was restricted 
to a head-on collision. Recently, Shibata and Uryu 
(2006) reported on first results of neutron star black 
hole binaries in full general relativity. 
Here, we will report mainly on our own, either New- 
tonian or pseudo- Newtonian results. These calcu- 
lations focused on the microphysics rather than the 
strong- field gravity aspect, therefore, the results con- 
cerning gravitational waves should be taken with a 
grain of salt. Details about the calculation of the 
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wave forms can be found in Rosswog et al. (2004), 
we will mention in the appropriate places in which 
direction fully relativistic calculations are expected 
to change the results. 

2.1. Double neutron stars 

For many years, neutron star masses, at least 
those in double neutron star systems, were thought 
to be tightly clustered around a mass of m= 1.35 
M (Thorsett and Chakrabarty 1999). Recent data, 
see e.g. Stairs (2004), show considerable deviations 
from this value for several of the known binary sys- 
tems. The recently discovered double neutron star 
system PSR J1756-2251 (Faulkner et al. 2005), for 
example, has a mass ratio of q= 0.84 (mi = 1.40 
and rri2 = 1.18 M ) and the double neutron star 
J1518+4904 even has q = 0.67 (mi = 1.56 and 
rri2 = 1.05 M ), but still relatively large errors (see 
Stairs 2004). 

For our new DNS simulations, we take masses of 
twice 1.4 M as the typical case, but we also con- 
sider the more extreme case with a 1.1 and a 1.6 
M neutron star. The dynamical evolution during 
the last moments of a binary system with twice 1.4 
M is shown in Fig. [H The neutron stars have no 
initial spin, their initial separation is just outside the 
plunging regime for our equation of state at 48 km 
and each star is modelled with slightly more than a 
million SPH particles (run DA in Tab.[T]). Once the 
plunging sets in, the stars merge within about one or- 
bital revolution. Excess angular momentum is shed 
into symmetric spiral arms that subsequently spread 
into an extended accretion disk, see Sec. 03 around 
a super-massive, neutron star-like object. The evo- 
lution of the 1.1 and 1.6 M system (run DB in 
Tab.[Q), similar to J1518+4904, is shown in Fig. O 
In this case the lighter star starts to transfer mass 
towards and then completely engulfs the more mas- 
sive component. The mass transfer leads to a com- 
plete tidal disruption of the 1.1 M star, parts of 
which form a tidal tail, see panel three of Fig. [2l 
For the chosen system parameters, the circulariza- 
tion radius of the accreted material is smaller than 
the accretor radius, therefore, like in an Algol sys- 
tem, the accretion stream directly impacts on the 
accretor crust, see panel one in Fig. [2] and Fig. 03 
For less extreme mass ratios the accreted material 
slides more smoothly around the surface of the heav- 
ier companion, otherwise the evolution is similar. 
In the last panels of Figs. [T] and [2] we show the re- 
tarded gravitational wave amplitudes, h+ and h Xl 
as seen along the binary axis by a distant observer 
at distance d. In both cases one sees the last stages 
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Fig. 4. Exponent of the mass-radius relationship for a 
polytropic star. 

of the "chirp" -signal up to about 2.5 ms, when the 
stars come into contact. For the chosen initial sep- 
aration of the 2 x 1.4 M case the peak frequency 
of « 1 KHz is reached 2.5 ms after the simulation 
start. If femin labels the minimum detectable ampli- 
tude, such a system would be visible out to a distance 
d = 18 Mpc (10 _21 //i m in). During the chirp phase 
up to about 2.5 ms the signal is the result of the 
binary orbital motion, after that, up to about 6 ms 
it is determined by the elongated central object and 
the spiral arms. The central object does not become 
rotationally symmetric by the end of the simulation 
and, therefore, keeps emitting quasi-periodic grav- 
itational waves at about a tenth of the amplitude. 
This is behaviour is expected for a stiff equation of 
state (EOS) such as the Shen et al. (1998) EOS 
that we use. It had been noted very early on (Ra- 
sio and Shapiro 1994) that this post-merger gravita- 
tional wave signal is a way to constrain the equation 
of state. 

For the 1.1-1.6 M system the peak frequency is 
slightly lower, ^ pea k ~ 780 Hz, the same statement 
is true for the peak amplitude. 

This basic picture of the dynamical evolution car- 
ries over to general relativistic calculations. The 
stronger gravity, however, produces neutron stars of 
higher compactness and a more compact remnant. 
Apart from the most obvious consequence, the pos- 
sible gravitational collapse to a black hole, a general 
relativistic merger yields higher peak frequencies and 
gravitational wave amplitudes and therefore ensures 
detect ability out to larger distances. For more de- 
tails we refer to the recent calculations of Shibata et 
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Fig. 5. Dynamical evolution of a neutron star black hole system. Left column: example of a system that undergoes 
episodic mass transfer (M ns = 1.4 Mq, Mbh= 3 Mq). A mini-neutron star keeps orbiting the hole, thereby transfers 
mass directly into the hole and sheds mass which subsequently spreads into a dilute disk surrounding the orbiting binary 
system. It is only after 47 orbital revolutions (t ~ 220 ms) that the mini-neutron star is finally disrupted, see panel 
4, and its remains form an accretion disk of ~ 0.05 M©, see Fig. 1111 This episodic mass transfer is imprinted on the 
gravitational wave signal, see Fig. [9] Right column: binary in which the neutron star is directly disrupted (M ns = 1.4 
Mq, Mbh= 14 M©). The figures in each column refer to 3, 5 and 37 ms, respectively, after the onset of mass transfer. 
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al. (2005). 
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Fig. 6. Dynamical evolution of a neutron star black hole 
system that undergoes episodic mass transfer; continua- 
tion of Fig. [5] left column. 



2.2. Neutron star black hole 



The accretion and merger dynamics of neutron 
star black hole systems is, mainly due to the larger 
possible mass ratios, substantially more complicated. 
It is a result of the interplay between gravitational 
wave emission, mass transfer and the nuclear EOS 
(Rosswog et al. 2004). Gravitational wave emission 
will drive the binary towards coalescence, the mass 
transfer from the lighter neutron star to the heavier 
black hole has the opposite tendency. The influence 
of the EOS is twofold. On the one hand it determines 
the radius of the star, i? ns , and therefore the tidal 
radius, at which the neutron star will be disrupted. 
A simple estimate for this radius is given by 



Rtid — 



Mbh 



1/3 



(i) 



On the other hand the EOS is responsible for the way 
the star reacts to mass loss, i.e. whether it expands 
or contracts. For a polytropic star the mass-radius 

r-2 

relationship is R oc M 3r - 4 , where T is the polytropic 
exponent of the EOS (e.g. Kippenhahn and Weigert 
1990). The exponent of the mass-radius relation of a 
polytropic star, (r - 2)/(3r - 4), is plotted in Fig.gJ 
For T < 2 the star expands on mass loss and for 
T > 2 it shrinks on mass loss. For a realistic EOS, 
however, the situation is complicated by the fact that 
T changes with density. For the EOS we use (Shen 
et al. 1998), the neutron star will shrink on mass 
loss for masses above 0.4 M and expand otherwise. 
To build up an accretion disk a large ratio of tidal 
radius to innermost stable circular orbit, Ri sco , 



R 



■tid 



R-n 



Rit 



6G (M ns M b 2 h )i/3' 



(2) 



is favorable. Thus, this simple estimate suggests that 
low masses for both the black hole and the neutron 
star are favorable. 

Taking together the calculations described in Ross- 
wog (2005) and those of Tabled] the whole plausible 
black hole mass range from 3 to 20 M (Fryer and 
Kalogera 2001) has been covered. For black holes be- 
low 10 M Newtonian potentials and for more mas- 
sive BHs Paczyhski -Wiita (PW) pseudo-potentials 
have been used. To mimic the presence of a last sta- 
ble orbit in the Newtonian calculations, we set up 
an absorbing boundary at Ri SCO = 6GMbh/c 2 . For 
the PW-potential matter is absorbed at 3GMbh/c 2 , 
for details see appendix in Rosswog (2005), but this 
material has passed the last stable orbit already and 
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falls freely towards the hole, so the exact location 
of the boundary is unimportant. In all calculations 
a non-spinning neutron star mass with 1.4 M was 
used. 

For illustration, we want discuss two of the runs in 
more detail: A run with a 3 M black hole (run NA 
in Tab. [1]) as an example for a long-lived, episodic 
mass transfer, and a 14 M run in which the neu- 
tron star is disrupted at during the first encounter. 

2.2.1. Episodic mass transfer 

In this calculation with a 3 M black hole (run 
NA in Tab. [T]) the neutron star undergoes a long- 
lived mass transfer without being disrupted, see 
Fig. left column. After a primary mass transfer 
episode with peak mass transfer rates beyond 100 
M /s, see Fig. [7J left panel, the neutron star re- 
cedes again and quenches the mass transfer. During 
this episode its mass is reduced to about 0.5 M , 
see right panel of Fig. [HI Subsequently, it under- 
goes several more mass transfer episodes, after each 
episode mass transfer practically shuts off, leaving 
the neutron star constant for a short time (see "step- 
like" behaviour in the right panel of Fig. [8]). After 
about 25 episodes the mass transfer settles into a 
stationary state in which the neutron star transfers 
approximately 1 M /s directly into the hole with- 
out forming an accretion disk. Its general tendency 
during this episode is to move further away from the 
black hole, see left panel of Fig. [HI Some material is 
tidally shed from the neutron star into a dilute disk 
around the still orbiting binary. It is only after 47 
orbital periods, at t~ 218 ms, that the orbiting mini- 
neutron star becomes tidally disrupted and forms an 
accretion disk around the hole, see panels two and 
three in Fig. [U This final tidal disruption sets in 
when the mini-neutron star reaches its lower mass 
limit of about 0.18 M , see right panel of Fig. [8] (for 
the mass radius relationship of the Shen et al. EOS 
see Fig. 1 in Rosswog et al. (2004)). 
This long-lived binary phase is imprinted on the 
gravitational wave signal, see Fig. [9l left panel. The 
signal "chirps" until mass transfer sets in. The 
primary mass transfer episode reduces the gravi- 
tational wave amplitude by more than a factor of 
two. While the neutron star is constantly stripped 
by the black hole, the gravitational wave amplitude 
decreases until the minimum neutron star mass is 
reached. During this phase the binary emits at a 
frequency of « 500 Hz, close to LIGO's sensitivity 
peak. This phase should be detectable out to a dis- 
tance of d = 10 Mpc (10 _21 //i m in). Once this stage 
is reached, the gravitational wave signal shuts off 



abruptly as the matter settles quickly into an ax- 
isymmetric disk. During this final disruption, the 
mass transfer rate increases for a last time beyond 
10 M /s, see the last pronounced peak around 200 
ms in the left panel of Fig. 

2.2.2. Direct disruption 

For more massive black holes, the self-gravity of 
the neutron star can be overcome already during 
the first, primary mass transfer. An example of 
such a case with a 14 M black hole is shown in 
Fig. right column. Although the neutron star 
is disrupted, a density peak remains visible in the 
debris spiral arm until the end of the simulation. In 
this case an accretion disk forms, but large parts 
of this disk are inside the innermost stable circular 
orbit of the black hole and therefore falling with 
large radial velocities into the hole. This disks are 
geometrically thin and, apart from a spiral shock 
from the self-interaction of the accreted matter, 
essentially cold. A detailed discussion of this 
calculation can be found in Rosswog (2005). During 
the first approach mass is transferred into the black 
hole at a peak rate of more than 1000 M /s, see 
right panel of Fig. 

Again the dynamics leaves a clear imprint on the 
gravitational wave signal: after the chirp stage the 
signal immediately dies off, see Fig. [9] right panel. 
At peak amplitude, the signal should be visible out 
to d = 250 Mpc (10- 21 /^min). 

For the higher mass black holes, we do not find 
accretion disks that are promising to produce ener- 
getic GRBs. While some systems do form low-mass 
disks, see Fig. [5j right column, for a 14 M BH, 
for black holes with M > 18 M , nearly the whole 
neutron star is directly fed into the hole, only a 
small fraction (< 0.08 M ) receives enough angular 
momentum to be ejected in form of a half-ring (see 
Figure 3 in Rosswog 2005) of neutron-rich debris 
material. 

2.3. Constraints on the nuclear equation of state 
from gravitational waves 

As mentioned earlier, the neutron star equation 
of state (EOS) could be seriously constrained by 
a detected gravitational wave signal. We want to 
briefly summarize gravitational wave signatures that 
carry information on the nuclear EOS. 

• A simple estimate for the neutron star radius 
can be obtained using Newtonian gravity and 
^gw = 2^ rb> where ^gw and ^ Q rb are the grav- 
itational wave and orbital frequencies. If one 
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assumes an equal mass binary, that the peak 
occurs once the neutron star surfaces touch and 
that changes due to tides are negligible, one 
finds 

/ gm y /3 

-„o, ( M V /lMHz\* , , 

= 16.8 km — — — .(3 

V2.8M ; \ v v ) K ' 

Here, u p is the peak GW- frequency at the end 
of the chirp phase and M the total mass of the 
binary. 

More accurately, using quasi-equilibrium se- 
quences of neutron star binaries, the compact- 
ness GM/Rc 2 can be determined from the final 
deviation of the gravitational wave energy spec- 
trum from a point mass binary signal (Faber 
et al. 2002). If additionally the neutron star 
masses are known from the inspiral signal (Cut- 
ler and Flanagan 1994), the neutron star radius 
can be derived and severe constraints on the 
equation of state can be imposed. 

• The stiffness of the EOS influences the post- 
merger behaviour of the remnant. If the EOS 
is stiff enough, the central object will keep a 
triaxial, ellipsoidal shape for many orbital pe- 
riods and thus will continue to emit gravita- 
tional waves (e.g. Rasio and Shapiro 1994), be- 
fore it (probably) finally collapses into a black 
hole. For a soft EOS it will quickly become axi- 
symmetric and the GW-signal will die off fast. 
The influence of quark matter on the gravita- 
tional wave signal has recently been investigated 
by Oechslin et al. (2004). 

• As was shown in Sec J2.2.1| NSBH systems with 
low black hole masses can undergo a long-lived 
(in comparison to orbital time), mostly episodic 
mass transfer activity. Test calculations with a 
soft (r = 2) polytrope (Rosswog et al. 2004) 
indicate that this behaviour is specific for stiff 
equations of state. Moreover, from a grav- 
itational wave signal such as the one shown 
Fig. [9j left panel, the time when the neutron 
star reaches its minimum mass can be directly 
read off: at this point the gravitational wave 
signal suddenly shuts off. 

3. ACCRETION DISKS 

The accretion disks produced in a compact bi- 
nary merger are in several respects different from 
a standard Shakura-Sunyaev (1973) accretion disk. 



Having just been built up in a violent disruption, 
they are far from a steady state and completely 
opaque to photons. Cooling is only possible via 
neutrino emission. The corresponding processes, 
however, are not generally fast enough to cool the 
disks on a timescale comparable to the dynami- 
cal timescales. The disks often span the full range 
from neutrino opaque to completely transparent with 
relatively large semi-transparent transition regions 
making an analytical treatment without restrictive 
assumptions difficult. Over the last years several 
groups have worked in different approximations on 
the neutrino-cooled disk problem (e.g. Popham et 
al. 1999, Narayan et al. 2001, Kohri and Mineshige 
2002, DiMatteo et al. 2002). Recently, equilibrium 
disks around rotating black holes have been con- 
structed by Chen and Beloborodov (2006). Neutrino 
opacity effects have been included in several of the 
recent numerical simulations (Ruffert et al. 1997, 
Rosswog and Liebendorfer 2003, Lee et al. 2005). 



3.1. Double neutron star mergers 

We show in Fig. [10] disk properties of run DA 
that we consider representative for the DNS case. 
The dominant neutrino-emitting parts of the disks 
have densities in the range from ~ 10 13 down to 
~ 10 10 gem -3 , see panel one, the temperatures 
in these regions are ~ 3 MeV. At the time of the 
merger, the neutron stars are still very close to cold 
/^-equilibrium as the tidal interaction is too short 
to cause substantial compositional changes at the 
prevailing temperatures (T < 10 8 K; Lai 1995). 
The weak interaction rates during the merger are 
sufficiently slow that the matter in the remnant is 
at ~ 15 ms still close to its initial electron fraction, 
Y e ~ 0.05, see panel three in Fig. [I0l Therefore the 
assumption of /^-equilibrium is not justified. 
The central object is apart from its surface layers 
completely opaque to the neutrinos. The inner parts 
of the debris torus still have large optical depths, 
see Fig. 8-10 in Rosswog and Liebendorfer (2003), 
only the outer parts are neutrino transparent. As 
the disks cannot cool on a dynamical timescale 
they are puffed-up and advection dominated. Inflow 
towards the central object proceeds mainly along 
the equatorial plane (see arrows in panel two and 
three) and along the disk surface, the flow inside 
the disks shows convective circulation. For a further 
discussion of the circulation patterns see Lee and 
Ramirez-Ruiz (2002), Rosswog and Davies (2002) 
and Lee et al. (2005). 



10 



ROSSWOG 




t [ms] t [ms] 

Fig. 7. Left panel: mass transfer (in M©/s) into the black hole for the 1.4 M© and 3 M© system with episodic mass 
transfer (run NA in Tab.[T|. Right panel: dito for the 1.4 Mq and 14 M© system in which the neutron star is directly 
disrupted. 




Fig. 8. Left panel: orbit of the mini- neutron star in the neutron star black hole system (M ns = 1-4 M© , Mbh= 3 M©) 
that undergoes episodic mass transfer. After mass transfer has set in, the neutron star orbits the black hole for 47 
orbital revolutions before it is finally disrupted. Right panel: elvolution of the neutron star mass (defined as the mass 
with p > 10 13 gem -3 ) of the same system. 
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Fig. 9. Left panel: gravitational wave signal of a binary system that undergoes episodic mass transfer (M ns = 1.4 M©, 
Mbh= 3 Mq). Right panel: gravitational waves for the 1.4 M© and 14 M© system in which the neutron star is directly 
disrupted (run II of Rosswog (2005)). 



3.2. Neutron star black hole mergers 

In the neutron star black hole case we find it 
much harder to build up massive accretion disks. In 
the low-mass cases the episodic mass transfer pre- 
vents disk formation until the neutron star is finally 
disrupted once it has been stripped down to its min- 
imum mass of about 0.18 M . Such a case was 
shown in detail in Figs. [5] and El The resulting disk 
as seen in panel six of Fig. [6] has about 0.05 M . For 
the cases with larger black hole mass, the neutron 
star is completely disrupted early on, see Sec. \2.2.2\ 
but close to the innermost stable circular orbit. As 
a consequence, matter falls with large radial veloc- 
ities into the hole within about one orbital period. 
These disks are essentially cold (apart from a spi- 
ral shock that occurs due to self-interaction of the 
accretion stream, see Fig. 6 in Rosswog 2005), low 
in density and, contrary to the neutron star merger 
disks, they are geometrically thin. For even higher 
black hole masses the neutron star is nearly swal- 
lowed completely without any disk formation. For a 
further discussion we refer to Rosswog (2005). 
These result seems to depend on the nuclear equa- 
tion of state. Janka et al. (1999) who use the 
softer Lattimer-Swesty EOS find disks that are more 
promising for GRBs. 



To quantify the amount of present debris material 
for both the DNS and the NSBH cases, we use the 
mass that is gravitationally bound and has a density 
p < 10 13 gem -3 . The evolution of this debris mass 
is shown in Fig. [TH In the DNS cases the debris 
mass is around 0.3 M , for the black hole cases it 
is about one order of magnitude smaller, at maxi- 
mum « 0.05 M . In run NA, which is described in 
the mass transfer is constantly 



detail in Sec. 12.2. 1[ 
driving neutron star oscillations, see Fig. [Til right 
panel. The debris mass only increases beyond 0.01 
M after the neutron star is completely disrupted at 
t« 218 ms. 

To further illustrate the different thermodynamic 
condition prevailing in the debris of DNS and NSBH 
remnants, we show in Fig. [12] trajectories in the 
p — T-plane. The left panel refers to the two neutron 
star cases, the right panel shows typical trajectories 
of NSBH systems. The solid lines always refer to 
the average of the hottest 10 % of the debris, the 
dashed lines show the average conditions of all mat- 
ter. While in the neutron star case the hottest mate- 
rial is always at very high densities (> 10 14 gem -3 ), 
the hottest regions in the NSBH cases move during 
the evolution down to densities of < 10 10 gcm _3 and 
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Fig. 10. Vertical structure of the debris disk resulting 
from an equal mass (1.4 M©) double neutron star merger 
(run DA, see Table 1, at t= 16.88 ms). Panel one shows 
the global mass density distribution in the XZ-plane; 
panel two zooms into inner disk region, color-coded is 
the temperature. To enhance the temperature contrasts 
in the disk the upper limit of the colorbar has been fixed 
to 6 MeV. Panel three shows the distribution of the elec- 
tron fraction, Y e ; in panels two and three the velocity 
field is overlaid. 

below. In the direct disruption case that was dis- 
cussed in Sec. I2.2.2[ the average temperature of even 
the hottest 10 % of the material quickly drops below 
1 MeV (solid green curve). 

4. NEUTRINO EMISSION 

4.1. Double neutron stars 

The debris in the DNS case is very neutron-rich 
and hot, so the neutrino luminosities are generally 
dominated by electron anti-neutrinos from positron 
captures onto free neutrons, e + + n — > p + v e . 
Electron neutrinos, mainly from electron captures, 
e + p — > n + v e , are second most important and 
followed by the heavy lepton neutrinos , T>^ , v T , v T 
that we generally refer to as v x . While the v e and 
v e are emitted predominantly from the inner disk re- 
gions or the surface layers of the central object, the 
heavy lepton neutrinos v x are mainly produced in 



the hot high-density regions inside the central object. 
Since the heavy lepton neutrinos are not absorbed 
by nucleons, the surrounding matter is more trans- 
parent to them. Therefore, they can escape from the 
hotter, high density regions more easily, which is why 
they have the largest average energy. The average 
energies of the neutrino species are relatively robust 
against changes in the system parameters, typically 
{E Ve ) « 7 MeV, (E Pe ) « 12 MeV and (E Ux ) « 23 
MeV. In the unlikely case that neutrinos from a neu- 
tron star merger would be detected, the substantial 
drop in the v x , but only a milder drop in the v e and 
v e luminosity would indicate the moment when the 
central objects collapses into a black hole. 
The neutrino luminosities after a DNS merger in- 
crease smoothly until a stationary state is reached, 
see Fig. [13l At this stage the luminosities of the 
numerical models can be roughly fit by 



Li 



M 



2.6 M P 



(4) 



with L\ = 6.8 • 10 52 erg/s, a = 4.5 and M being the 
total binary mass. This formula does not apply to 
cases where an Algol-like impact occurs (run DB), in 
such a case the neutrino-luminosities are higher due 
to the very strong shock heating at moderate densi- 
ties. The contribution of the different neutrino fla- 
vors to the total luminosity changes with the details 
of the system under consideration, but they are ap- 
proximately given by L Ve « 0.3 • L* ot , L Ve « 0.5 • L l ° l 
andL^ «0.2.L£ ot . 

The disks that were produced in neutron star black 
hole calculations were only of moderate masses and 
temperatures and therefore produced much lower lu- 
minosities than in the neutron star case. Calcula- 
tions with the softer EOS of Lattimer and Swesty 
(1991) (Janka et al. 1999), find results not too dif- 
ferent from the DNS case. 

5. GRBS 

Compact binary mergers have been recognized 
as possible central engines of GRBs since many 
years (Blinnikov et al. 1984, Eichler et al. 1998, 
Paczyriski 1991, Narayan et al. 1992). For general 
reviews on the manifold aspects of gamma-ray 
bursts we refer to recent reviews (e.g. Piran 2005 or 
Meszaros 2006). 

To attain the ultra-relativistic motion required to 
explain both the short-time variability and the 
non-thermal spectrum of GRBs, large amounts of 
energy have to be deposited in a region of space 
that is practically devoid of baryons. Two popular 
mechanisms to achieve this are the annihilation of 



LAST MOMENTS 



13 




Fig. 11. Comparison between double neutron star and neutron star black hole mergers. Left panel: debris mass 
(gravitationally bound, p < 10 13 gem -3 ) for the double neutron star case. Right panel: debris mass for the neutron star 
black hole cases. Note the different scales on the axes of both panels. 




Fig. 12. Comparison between double neutron star and neutron star black hole mergers. Left panel: Trajectories in the 
p — T-plane for double neutron star cases (run DA and DB), right panel: typical neutron star black hole cases. For each 
case, the trajectories of the hottest 10% of the material and an average over all matter are shown. 
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Fig. 13. Neutrino luminosities of the double neutron star 
merger cases. Solid lines refer to the case with 2 x 1.4 
Mq, the dashed lines to the 1.1 and 1.6 M0case. 

neutrino anti- neutrino pairs (Eichler et al. 1989, 
Mochkovitch et al. 1993, Ruffert et al. 1997, 
Popham et al. 1999, Asano and Fukuyama 2003, 
Rosswog and Ramirez-Ruiz 2002, Rosswog et al. 
2003) and mechanisms that rely on (ultra-)strong 
magnetic fields (Narayan et al. 1992, Usov 1992, 
Duncan and Thompson 1992, Thompson 1994, 
Meszaros and Rees 1997, Kluzniak and Ruderman 
1998, Lyutikov et al. 2003, Rosswog et al. 2003). 

5.1. Vi — annihilation 

The energy deposition rate per volume via neu- 
trino anti-neutrino annihilation (see Ruffert et al. 
1997) can be written is discretised form as 



i=e,fi,r 

= E ^E^eIm^ + ^VL 



z=e,fi,T 



L k 



L k [ {E Vi ) k + (E^f 
4 if 4, {E Vi ) k {E 9i ) 



where the index i labels the type of neutrino. Here 
L k is the neutrino luminosity of grid cell fc, is the 
distance from the centre of grid cell k to the point r, 
dk = \r— r)fe|, {E Vi ) k is the average neutrino energy in 
grid cell /c, Hkk' — 1 — cos 6kk' and Okk' is the angle at 
which neutrinos from cell k encounter anti-neutrinos 
from cell k' at the point r. The constants are given 



T^T^Y^v ~ C A f + (Cy + C A )% 

A h » = A hr = ^^^[(Cy - C A f + (Cy + 

C A - 2) 2 ], A 2 , e = ^^[2C 2 V - C% A 2 ^ = A 2 , r = 
^^[2(C V - l) 2 - (C A - l) 2 ], where Cy = 1/2 + 
2 sin 2 W , C A = 1/2, sin 2 Q w = 0.23 and a = 1.76 x 



10- 44 cm 2 . 



The thick-disk geometry that is a natural out- 
come of a double neutron star merger, see Fig. [101 
is favorable for an efficient annihilation as neutrinos 
have enhanced probability to collide head-on. The 
most interesting regions for the energy deposition 
are the centrifugally baryon-cleaned regions along 
the original binary rotation axis. Contours of 
the matter density and the annihilation energy 
deposition rate can be found in Fig. 2 in Rosswog 
et al. (2003). The deposition of a large amount of 
energy in a photon-pair plasma between the funnel 
walls results in a bipolar outflow along the binary 
rotation axis. The ratio of deposited annihilation 
and rest mass energy determines the asymptotic 
Lorentz-factors of the resulting outflow. For a 
double neutron star merger these Lorentz-factors 
lie between a few dozens and a few times 10 4 (see 
Rosswog et al. 2003, Figs. 4-6). The total energy 
contained in this relativistic outflow is ~ 10 48 erg 
(Rosswog and Ramirez-Ruiz 2003). 
Similar to the case of a new-born proto-neutron 
star, the huge neutrino luminosity will also drive 
an energy baryonic winds off the remnant. The 
interaction with the baryonic material collimates 
the neutrino-annihilation-driven bipolar outflow. 
Using simple estimates Rosswog and Ramirez-Ruiz 
(2003) found a broad distributions of opening angles 
and isotropised luminosities, centered around 7° 
and a few times 10 50 erg, respectively. 
More recently, Aloy et al. (2005) have performed 
detailed, axisymmetric relativistic hydrodynamic 
simulations in which they parametrised the energy 
deposition above remnant accretion disks. Some of 
their models did not produce ultra-relativistic out- 
flows and instead lead to low-luminosity UV-flashes, 
others were able to produce ultra-relativistic, struc- 
tured outflows that are promising GRBs sources. 
While the ultimate outcome seems to depend quite 
sensitively on the details of the considered system 
(disk structure, energy deposition rates, injection 
duration), their most promising systems produced 
ultra-relativistic outflows with half-opening angles 
of 5-10° and apparent isotropic energies of up to 
10 51 erg. 

It has to be stressed that to date no self-consistent 
calculations of a full merger with neutrino emission, 
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Fig. 14. Neutrino luminosities of the neutron star black hole cases. Left panel: black holes with 3 (dashed) and 4 
Mq (solid curves). Right panel: black holes with 7 (solid) and 10 Mq (dashed curves). 



annihilation, launch of winds and relativistic out- 
flow could be performed. Only individual aspects 
have been simulated that are a posteriori patched 
together into a global picture. Nevertheless, the 
theoretical numbers found in the last two investi- 
gations are consistent with recent observations of 
short GRBs that find values of E 1 of a few times 
10 48 erg and opening half-angles from about 9 to 
about 14° (Fox et al. 2005, Berger et al. 2005). 
Note however, that there are recent observations 
(Berger et al. 2006b) that suggest that a fraction of 
short GRBs may have substantially larger energies 
than the first few short bursts that were detected. 
Thus, neutrino annihilation from compact merger 
remnants, at least in the DNS case, seems to be in 
reasonable agreement with recent observations of 
short GRBs. 

5.2. Magnetic fields 

Nevertheless, it is hard to see how magnetic 
fields should not play an important role. Most neu- 
tron stars are endowed with strong magnetic fields: 
young pulsars typically have ~ 10 12 G and magne- 
tars (Duncan and Thompson 1992) are thought to 
have magnetic fields between 10 14 and 10 15 G. While 
still more speculative than the neutrino annihilation 
mechanism in the sense that the existing calculations 
are less detailed, there are good arguments to expect 
a dramatic growth of existing seed fields. 
The merger debris is differentially rotating, both the 
central object before collapse, see Fig. 9 in Rosswog 
and Davies (2002), and the surrounding disk. Thus 
the initial seed fields can grow either just by a lin- 
ear winding up the field lines or, more likely, via 



the magnetorotational instability (MRI) (e.g. Bal- 
bus and Hawley 1998). Recent magnetohydrody- 
namics simulations (Price and Rosswog 2006) sug- 
gest that magnetic seed fields are rapidly amplified in 
the Kelvin-Helmholtz instability at the shear inter- 
face between the two neutron stars, see Fig. [151 This 
amplification occurs on a timescale of only ~ 1 ms, 
long before the collapse can set in. As the shortest 
modes grow fastest in a Kelvin-Helmholtz instabil- 
ity, the reached maximum field strengths are limited 
by the finite numerical resolution. At the highest af- 
fordable resolution, about 2 • 10 15 G are reached, but 
very plausibly much higher field strengths will be re- 
alized in nature. The kinetic energy of the central 
object is large, about 8 • 10 52 erg, and if say 10 % of 
this energy can be transformed into magnetic field 
energy, the field strength averaged over the central 
object will be in excess of 10 17 G. Locally, the field 
could be substantially higher. 

To become dynamically important, the fields have 
to become extremely strong. The field strengths at 
which the magnetic pressure balances the nuclear 
matter pressure are shown in Fig. [TU In the central 
object of a merger these field strengths are in excess 
of 10 1T G, in the accretion disk they are still be- 
tween 10 14 and 10 16 G. If the central object remains 
stable for long enough, it could provide a Poynting- 
dominated, scaled-up relativistic pulsar wind. In this 
case energy is released at a rate of 



(f ) md ~10 51 J B 1 2 6 P™X 6 erg/ S (6) 
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Fig. 15. Velocity fields during the merger of two 1.4 
Mq neutron stars (run DA). The shear interface, see 
panel one, becomes Kelvin-Helmholtz unstable and forms 
a string of vortex rolls (panel two, velocities in corotating 
frame) in which the magnetic field is curled up. 

and the central object will spin down on a timescale 
given by 



Tsd 



zs d 16 r- ms 



20 km 

Rco 



2.5 M 



(?) 



where B 16 is £/10 16 G and P ms = P/l ms. 

Fluid instabilities, such as the Kelvin-Helmholtz 
instability illustrated for run DA in Fig. [151 will lo- 
cally curl up the magnetic field lines. This had been 
suggested earlier (Rosswog et al. 2003) and has re- 
cently been confirmed by first MHD-merger simula- 




Fig. 16. Field strength (in G), where the magnetic pres- 
sure, equals the matter pressure (calculated from the 
Shen et al. EOS (1998) with T = 0.5 MeV and Y e = 0.1). 



tions (Price and Rosswog 2006). Once such high- 
field pockets reach field strengths close the equipar- 
tition, see Fig. [161 they will become buoyant, float 
up and produce explosive reconnection events (e.g. 
Kluzniak and Ruderman 1998, Rosswog et al. 2003). 
This may also occur in the disk (e.g. Narayan et al. 
1992), though at a lower field strength, see Fig. [T6l 
Once the central object has collapsed, the "standard 
GRB-engine", a black hole plus accretion torus, re- 
mains, for which various energy extraction mecha- 
nisms have been suggested (e.g. Blandford and Zna- 
jek 1977, Blandford and Payne 1982 or McKinney 
2005). For a discussion of these mechanisms we re- 
fer to the literature. 

As discussed in the previous sections in our neutron 
star black hole calculations we find it difficult to form 
massive disks. The neutrino luminosities are even in 
the more optimistic cases two orders of magnitude 
lower than in the DNS case, see Figs. [13] and [TH 
This is mainly due to the much smaller debris mass, 
see Fig. [Til As the neutrino annihilation roughly 
scales with the square of the neutrino luminosity, see 
Eq. ([5]), neutrino annihilation deposits substantially 
less energy than in the DNS case. NSBH binaries 
would therefore not produce energetic GRBs, but 
rather a low-luminosity tail to the short GRB distri- 
bution. Possible observable signals are discussed in 
Rosswog (2005). 

We want to stress again that these results are sensi- 
tive to the nuclear EOS. With the softer Lattimer- 
Swesty-EOS Janka et al. (1999) find more promising 
conditions. 
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6. FORMATION OF A MAGNETAR? 

There are essentially two proposed formation 
mechanisms for pulsars: i) the collapse of very 
rapidly rotating stars with ordinary magnetic fields 
or ii) the collapse of extraordinarily magnetized stars 
with ordinary rotational speeds. In the first case 
the new-born proto-neutron star has to be rotating 
rapidly enough so that rotational periods are smaller 
than the overturn times, ~ 3 ms, of convective eddies 
(Duncan and Thompson 1992, Thompson and Dun- 
can 1993, Thompson and Duncan 1995), so that an 
efficient a-c<;-dynamo can become active. In this way, 
fields as strong as ~ 10 15 G can be plausibly gen- 
erated. Such a new-born neutron star will be spun 
down quickly by magnetic dipole radiation (and pos- 
sibly mass loss) and this will fuel the supernova ex- 
plosion with an additional 2tt 2 I/P 2 - 2-10 52 / 4 5/^m S 
ergs stemming from rotational energy. Here we have 
used I45, the moment of inertia in units of 10 45 gcm 2 
and P ms , the rotation period as measured in millisec- 
onds. Therefore, particularly bright supernova ex- 
plosions going along with magnetar formation are a 
natural prediction of this formation path. However, 
the three supernova remnants that seem to be associ- 
ated with magnet ars are remarkably unremarkable: 
rather than being particularly energetic they only 
have kinetic energies of ~ 10 51 ergs (Sasaki et al. 
2004, Vink 2006). An alternative formation scenario 
could be that the magnetar field strengths are ob- 
tained via magnetic flux conservation alone during 
the collapse of highly magnetized progenitor stars 
(Ferrario and Wickramasinghe 2005). 
Here, we want to discuss a third possible formation 
channel, namely a neutron star coalescence that pro- 
duces produces a stable, magnetar-like remnant. 

The first point to clarify is whether the DNS 
merger rates are large enough to possibly contribute 
a non-negligible fraction of magnetars. The esti- 
mated rates at which double neutron systems (DNS) 
merge have been constantly increasing over the last 
years. The discovery of the double pulsar PSR 
J0737-3039A (Burgay et al. 2003) alone has in- 
creased the estimated merger rate based on obser- 
vational data by an order of magnitude. Current 
estimates (Kalogera et al. 2004) are in the range 
from 4 to 224 -10 -6 per year and Galaxy. Nakar 
et al. (2005) based their estimates on the assump- 
tion that short-hard gamma-ray bursts result from 
DNS mergers. They find a "best guess" value that 
is higher than the above upper limit by more than 
an order of magnitude. This last rate may be some- 
what optimistic, but it is fair to state that the DNS 
merger rate is comparable to the magnetar birthrate 



in our Galaxy, about 1-10 -10 -4 per year (Duncan 
and Thompson 1992). DNS mergers could therefore 
plausibly contribute to the magnetar production. 
The main question that needs to be addressed is 
whether the final remnant of a DNS merger can have 
a mass that can be permanently sustained by the 
equation of state. The maximum neutron star mass 
is unfortunately still only poorly constrained. An 
upper limit based on very general physical princi- 
ples such as causality is 3.2 M (Rhoades and RufEni 
1974)0. In a review of realistic models of nuclear 
forces Akmal et al. (1998) find that realistic models 
of nuclear forces limit the maximum mass of neutron 
stars to be below 2.5 M . The masses with the small- 
est error bars come from the observation of DNS and 
they are consistent with a remarkably narrow un- 
derlying Gaussian mass distribution with M= 1.35 
± 0.04 M (Thorsett and Chakrabarty 1999). Ob- 
servations of neutron stars in binary systems with a 
white dwarf, however, yield consistently higher neu- 
tron star masses, but the error bars are to date still 
much larger than in DNS. The most extreme case 
is PSR J0751+1807 with an estimated pulsar mass 
of 2.1 ± 0.2 M (Nice et al. 2005). Recent obser- 
vations of the neutron star EXO 0748-676 also seem 
to support large possible neutron star masses and 
have been interpreted as hints that neutron stars 
may have a conventional neutron-proton composi- 
tion (Ozel 2006). Here, we want to explore the pos- 
sibility that the final outcome of a DNS coalescence 
is, after some mass loss, a very massive and highly 
magnetised, but stable neutron star. 
Let us first estimate the gravitational mass of a 
merger remnant. Lattimer and Yahil (1989) found 
empirically that 



with a = 0.084 M , yields a good fit of the binding 
energy, B, as a function of the gravitational mass 
M. Thus, a standard binary system of twice 1.35 
M corresponds to a baryonic mass of about 3.0 M , 
or, a merger remnant made of all the baryons of the 
initial binary system would have 2.48 M of gravi- 
tational mass. 

We plot in Figure \17\ the baryonic mass that needs 
to become unbound in order for the final remnant 
to be above a specified maximum neutron star mass 
M max : 

m loss = 2(«M n 2 s + M ns ) - («M^ ax + M max ), (9) 

3 See Psaltis (2004) for a critical analysis of the underlying 
assumptions. 
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where we have restricted ourselves to the case of 
equal mass neutron stars. For example, if we assume 
the maximum neutron star mass to be 2.2 M , so 
just beyond the estimated mass of PSR J0751+1807, 
a system with two 1.2 (1.4) M neutron stars would 
need to loose a baryonic mass of 0.035 (0.52) M . 
Mass can be lost by dynamical ejection or by the 
combined action of neutrinos, rotation and strong 
magnetic field as a neutrino-magnetocentrifugally 
driven outflow. The latter mechanism has been dis- 
cussed recently (Thompson et al. 2004) in the con- 
text of the formation of magnet ars in supernovae. 
The amount of mass that is dynamically ejected is 
sensitive to the nuclear equation of state (EOS) with 
stiffer EOSs ejecting more than softer ones (Ross- 
wog et al. 2000). With the EOS of Shen et al. 
(1998), which is at the stiffer end of the EOS spec- 
trum and therefore consistent with the assumption 
of a relatively large maximum neutron star mass, we 
find that typically m e j ^0.03 M are dynamically 
ejected. 

The combined action of neutrino-heating, strong 
magnetic field and rapid rotation is hard to quantify 
as details of the (unknown) magnetic field geome- 
try enter and the mass loss rates depend sensitively 
on the neutrino luminosity and the rotational pe- 
riod (Duncan et al. 1986, Qian and Woosley 1996; 
Thompson et al. 2001 and Thompson et al. 2004). 
Strong magnetic fields enforce near-corotation of the 
wind with the stellar surface out to the Alfven ra- 
dius, Ra, where the kinetic enery, pv 2 /2 approxi- 
mately equals the magnetic energy density B 2 /Stt. 
Thus, wind mass elements are, like spokes on a bike 
wheel, forced to corotate with the remnant, lead- 
ing much faster loss of angular momentum and spin- 
down. Thompson et al. (2004) find that a strong 
magnetic field together with the rapid rotation will 
drastically increase the mass loss rate, in the fast- 
rotation limit they find a functional dependence of 



M oc exp(cj 2 ), 



(10) 



where uo is the angular frequency. This indicates that 
a substantial mass loss is possible. For example, for 
a central remnant period of 0.8 ms (see Figure 10 
in Rosswog and Davies 2002) and a neutrino lumi- 
nosity of Lp e = 8 • 10 51 erg/s they find a mass loss 
rate M ~ 0.1 M /s. If we assume that the results 
scale like M ex L^ 5 M~ 2 (Qian and Woosley 1996) 
and use the empirical relation Eq. (|4]), we can esti- 
mate M as a function of the total mass of the DNS, 
M. This is shown as the solid black line in Fig. [TT] 
for a duration of r = 0.5s, the total unbound mass 
being m ~ M • r + m e j . Taking these numbers at 




Fig. 17. Shown is the baryonic mass that has to be lost 
from a double neutron star system with twice M ns if 
the maximum gravitational neutron star mass is A^fmax • 
All masses are in solar units. The black solid curve is 
the estimate of the mass that is lost by the neutrino- 
magnetocentrifugal wind plus dynamic ejection, see text. 



face value would mean that, if the maximum neutron 
star mass is M max = 2.2 M (violet line in FigJTYj) , 
then even a standard DNS with twice 1.35 M would 
loose enough mass to leave a stable remnant. For 
M max = 2.0 M (blue), neutron stars with or below 
1.15 M would still produce a stable neutron star. 
To summarize, although hard numbers are difficult 
to obtain, the lower end of the DNS mass distribu- 
tion could plausibly produce stable, highly magne- 
tized neutron stars rather than black holes. 



7. LATE-TIME CENTRAL ENGINE ACTIVITY 

It was one of the surprises of the recent after- 
glow detections of short GRBs that GRB 050709 and 
GRB 050724 exhibited long-lasting X-ray flares that 
occurred at least 100 s after triggering the bursts. 
While the observations have been claimed to be sug- 
gestive for compact mergers, it is not immediately 
clear whether the merger model can accommodate 
this late time X-ray flaring activity. The standard 
central engine, black hole plus disk, has -at least in 
its most simple form- difficulties to accommodate 
these long time scales. The dynamical timescale for 
an accretion disk is 



Td 



0.05s 



K 



R 



200 km 



(U 

a 



2.5M, 



M bh 



(11) 

where a is the Shakura-Sunyaev viscosity parameter 
(Shakura and Sunyaev 1973) and Qk the angular 
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frequency of a ring of matter at distance R. Ob- 
viously, the estimated accretion timescale is much 
shorter than the timescales at which the X-ray 
flaring activity is observed. 

The X-ray activity might be caused by either a 
long-lived central engine, or by a wide distribution 
of Lorentz factors, or else by the deceleration of a 
Poynting flux dominated flow. Various possibilities 
are discussed in detail in Zhang et al. (2005). Here 
we want to address two possibilities: i) that the 
central object from a DNS merger survives for long 
enough to produce the flaring and ii) we also address 
the fallback accretion luminosity in the aftermath 
of a compact binary merger. 



7.1. Survival of the central object in a double 
neutron star merger 

The continued central engine activity could be 
due to the survival of the central object, either in a 
meta-stable state or -as discussed above- the rem- 
nant may in some cases loose enough mass to pro- 
duce a stable, highly magnetised neutron star. 
A meta-stable state is likely to occur, stabilisation 
could come for example from differential rotation 
(e.g. Morrison et al. 2004) or trapped neutrinos to- 
gether with non-leptonic negative charges (Prakash 
et al. 1995), but the time scale until collapse is un- 
certain. This is because all the imprecisely known 
aspects of the problem such as the EOS at supra- 
nuclear densities, the remnant radius and magnetic 
field, ellipticity of the central object and so on, enter 
with large powers in the estimate of the time scale. 
Therefore in the meta-stable case a large range of 
time scales until collapse can be expected for differ- 
ent initial binary systems. With reasonable param- 
eters even time scales as long as weeks are possible 
(see Rosswog and Ramirez-Ruiz 2002). 
If the GRB is indeed caused by the combined ac- 
tion of neutrino annihilation and buoyant pockets 
of ultra-strong magnetic fields, this will plausibly 
also cause late-time flaring activity, provided that 
the remnant is stable for long enough. The mag- 
netic field amplification acts on a time scale of about 
a milli-second (Price and Rosswog 2006) and thus 
large magnetic fields can be built up before the neu- 
trino luminosity has reached its peak at « 20 ms 
after the merger, see Fig. [131 In the early stages a 
buoyant high-field bubble can expand into a practi- 
cally baryon-free environment and thus reach large 
Lorentz-factors. At neutrino luminosities in excess 
of 10 52 erg/s the remnant also drives a very ener- 
getic baryonic wind. If the right conditions are met, 



the baryonic material can help to collimate the out- 
flow (Aloy et al. 2005), but it also poses a potential 
threat to the emergence of an ultra-relativistic out- 
flow. Initially, as the rapidly spinning remnant pro- 
duces high-field pockets in rapid succession, mag- 
netic pressure may help to keep the funnel above 
the central object baryon-clean. But as the remnant 
is constantly braked by magnetic fields and gravita- 
tional wave emission, it takes longer and longer to 
reach buoyancy field strength, while more and more 
baryons are ablated from the remnant by neutrinos. 
The first explosive reconnection events escape into 
a very clean environment, those produced later will 
be released into an envelope of neutrino-magneto- 
centrifugally expelled baryons. This different baryon 
loading may be responsible for the difference between 
bursts and X-ray flares. 

7.2. Fallback accretion 

Some of the matter expelled by gravitational 
torques is still gravitationally bound to the central 
object, and will fall back towards the remnant (see 
for example Figure El panel three). Direct numerical 
calculations over the interesting timescales are not 
feasible, therefore we follow a simplified analytical 
approach. We assume that the motion of this ma- 
terial can be approximated as Keplerian motion in 
a central potential produced by the remnant. From 
the total energy, E{, and the angular momentum, J^, 
of each SPH-particle the eccentricity of the particle 
orbit, ei, can be calculated 



= 1 



2EjJ[ 
G 2 mfM 2 ' 



(12) 



Here, rrii is the particle and M the enclosed mass. 
For the particles with e$ < 1, the semi- major axes 
are calculated, a$ = —GmiM/(2E i ) 1 from which the 
distances of the closest and farthest approach follow: 
#min = a»(l - Ci) and i? max = a*(l + e*). A particle 
with a velocity V{ that is currently located at radius 
Ti will reach the radius R^is after a time of 



where I ri ,r 2 is given by 



for Vi • Ti > 
for V{ • ?i < 



(13) 
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with D = 4AC - B 2 and A = B = 2GM and 
C = -$ (Rosswog 2006). 

i 

For the radius, where the fallback energy is dissi- 
pated, i^dis, we take in the double neutron star case 
the disk radius at the end of the numerical simula- 
tion. This is conservative in the sense that the disk 
will shrink, see below, and therefore the assumption 
of a radius fixed at Rai S yields shorter timescales 
and lower accretion luminosities. For the black hole 
cases, I choose Rais = 10GM/c 2 , so just outside 
the innermost stable circular orbit of a non-rotating 
(Schwarzschild-) black hole at Risco — 6GM/c 2 . 
The details of the fallback times and energies change 
slightly with i?dis> but none of the conclusions de- 
pends on the exact numerical value of i^dis- 
Fig. [18] shows the accretion luminosities, L acc = 
dEfb/dt, derived for various DNS and NSBH sys- 
tems, for details see Rosswog (2006). Here, E^ de- 
notes the difference between the potential plus ki- 
netic energy at the start radius, n, and the po- 
tential energy at the dissipation radius, i^dis- The 
curves have been obtained by binning the energies 
contained in the fallback material, Efb, according to 
the corresponding fallback times, r^, see Eq. (fT3j) . 
A fraction e of this energy is channeled into X-rays, 
Lx = eL acc . 

The double neutron star cases are rather homoge- 
neous with respect to their fallback accretion, in 
all cases the fallback material is approximately 0.03 
After an initial, short-lived plateau the lu- 
minosity falls off with time close to the expected 
5/3-power law (|Reeslll988t IPhinnevI 1 1989h . The last 
point in these curves is determined by the numer- 
ical mass resolution in the hydrodynamics simula- 
tions and should therefore be interpreted with some 
caution. All other points should be a fair represen- 
tation of the overall fallback activity. Typically, the 
X-ray luminosity about one hour after the coales- 
cence is Lx ~ (or) ' 10 44 erg/s. For the investigated 
mass range the spread in the luminosities one hour 
after the coalescence is about one order of magni- 
tude. 

The neutron star black hole cases show a larger di- 
versity. The mass in the fallback material of different 
mass ratios varies by about a factor of 500, see Ross- 
wog (2006), an hour after the merger the accretion 
luminosities of the different NSBH systems differ by 
about two orders of magnitude. The involved time 
scales change strongly with the binary mass ratio. 
For example, the 1.4 and 4 M NSBH case does not 
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Fig. 18. Fallback accretion luminosity for various com- 
pact binaries, taken from Rosswog (2006). Circles re- 
fer to double neutron star systems (DNS), triangles to 
neutron star black hole mergers (NSBH). Note that the 
rightmost point is determined by the mass resolution of 
the simulation. For reference a straight line with slope 
5/3 is shown. 

produce much eccentric fallback material. Accretion, 
at least to the resolvable level, is over in ~ 0.2 s. 
This accretion period may produce a short GRB, 
but probably not much X-ray activity. The 1.4 and 
18 M NSBH system is at the other extreme: its 
peak luminosity is lower by three orders orders of 
magnitude but extends (at a resolvable level) up to 
about one hour. The mass ratios in-between could 
possibly produce a (weak) GRB and extended X-ray 
activity up to about one day after the burst. 

8. SUMMARY 

We have discussed various aspects of the last 
stages in the life of a compact binary. In particu- 
lar we have addressed 

• Gravitational waves as a probe of the merger 
dynamics 

We have shown how the merger dynamics is im- 
printed on the gravitational wave signal. Ex- 
amples include the peak gravitational wave fre- 
quency of a double neutron star merger and 
the subsequent ringdown emission which both 
are sensitive to the equation of state. For a 
stiff EOS, low-mass black hole neutron star sys- 
tems undergo a long-lived, mostly episodic mass 
transfer phase. We show an example where the 
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mass transfer continues for as long as 47 orbital 
periods. This last binary phase is responsible for 
a gravitational wave signal with slowly decaying 
amplitude and increasing frequency. The signal 
only shuts off when the neutron star is finally 
disrupted upon reaching its minimum mass. 
The coincident of detection of a short GRB to- 
gether with a "chirping" gravitational wave sig- 
nal would be the ultimate proof of a compact 
binary central engine. 

• Accretion disks 

Compact binary mergers produce in some (but 
not all) cases accretion disks from the neutron 
star debris. These disks differ from the stan- 
dard thin disk accretion model in various ways. 
Due to insufficient cooling by neutrinos they are 
often thick and puffed up, they are usually nei- 
ther in a dynamical nor in /^-equilibrium. For 
double neutron star mergers they are generally 
massive, of order 0.25 M©, with temperatures 
of several MeV and an electron fraction of close 
to 0.1 (at least initially). For neutron star black 
hole systems we find it substantially more diffi- 
cult to form massive disks. For low-mass black 
hole systems, the initial episodic mass transfer 
with a surviving neutron star prevents the build 
up of a massive disk. Only after a complete 
disruption of the neutron star a disk of « 0.05 
M forms. A comparison of the debris masses in 
double neutron star and neutron star black hole 
mergers is shown in Fig.[TTl We find the neutron 
star black hole systems to produce disks masses 
that are lower by at least a factor six. For higher 
mass black holes the neutron star is disrupted 
early on, but very close to the innermost sta- 
ble circular orbit. The resulting accretion disks 
are geometrically thin and essentially cold. For 
black hole masses larger that about 16 M no 
disk forms at all, most of the neutron star is fed 
directly into the hole, the rest is dynamically 
ejected. 

These results are sensitive to the equation of 
state, for a softer equation of state the neutron 
star is disrupted more easily. 

• Gamma-ray bursts 

Neutrino annihilation from the remnant of a 
double neutron star merger can plausibly pro- 
vide the driving stresses to launch highly rela- 
tivists, bipolar outflows, but a large diversity 
in the observable burst properties and moder- 
ate burst energies of ~ 10 48 are expected. It is 



conceivable that a good fraction of systems fails 
to provide the right conditions and instead pro- 
duces X-ray or UV-flashes. 
MHD simulations of double neutron star merg- 
ers show that the magnetic field grows within 
the first millisecond to (probably much) beyond 
magnetar strength, possibly allowing for a burst 
production a la Kluzniak and Ruderman (1998). 
The burst is most likely a result of the combined 
action of both neutrino annihilation and ultra- 
strong magnetic field. 

Of the neutron star black hole binaries only sys- 
tems with low black hole masses could possibly 
form disks that are hot and dense enough to 
launch GRBs (unless the true nuclear equation 
of state is much softer than the one we use here). 
Among them, rapidly spinning holes are pre- 
ferred as both the last stable orbit and the event 
horizon move closer to the hole. The high-mass 
end of neutron star black hole binaries proba- 
bly produces a low-luminosity tail of the short 
GRB-distribution. 

• Magnetar formation 

If the maximum neutron star mass should in- 
deed be as high as indicated by recent obser- 
vations the amount of mass that would need 
to be lost after a merger to evade the fi- 
nal collapse to a black hole is not implausi- 
bly large. Therefore we want to point out 
the possibility that binary systems from the 
lower end of the mass distribution could pro- 
duce a highly magnetised, but stable neutron 
star rather than a black hole. Magnet ars outside 
normal supernova remnants would be obvious 
candidates for this formation channel. If true, 
this channel should also produce large amounts 
of neutron-rich nuclei synthesised in both the 
dynamically ejected neutron star material and 
in the neutrino-magnetocentrifugal wind that 
made the survival of the remnant possible in 
the first place. Being synthesised in a neutrino 
bath dominated by electron anti-neutrinos the 
wind material will to be more proton-rich than 
the dynamically ejected material. The exact Y e 
is determined by both the expansion time scale 
and the ratio of neutrino and anti-neutrino lu- 
minosities, but is roughly (Qian and Woosley 
1996) 



Y e ~ 0.28 



1 



L P J5 • 10 52 erg 
L, e /2-10 52 erg 



(15) 



where L v and L„ e are the electron neutrino and 
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anti- neutrino luminosities, respectively. The 
more neutron-rich dynamical ejecta with Y e « 
0.1 would be expected to be located at larger 
distances from the magnetar than the wind ma- 
terial. 

• Late time X-ray activity 

In the speculative case that the central object 
of a neutron star merger survives, late-time flar- 
ing would be caused by the same mechanism 
that produced the burst. If the burst is (at 
least in part) produced by explosive reconnec- 
tion events of buoyant high-field bubbles as sug- 
gested by Kluzniak and Ruderman (1998), parts 
of the remnant that become buoyant at later 
times have to expand into an environment that 
is already heavily baryon-polluted by neutrino- 
driven winds. Due to the higher baryon loading 
the emission would be mainly in the X-ray band 
of the spectrum. 

Much less speculative is the accretion of fall- 
back material which occurs naturally in the af- 
termath of a compact binary merger. A few 
percent of a solar mass are ejected into eccen- 
tric, but still bound orbits. The fallback of this 
material produces extended X-ray activity up to 
many hours after the main burst. 

Some of the figures have been produced using the 
software SPLASH kindly provided by Daniel Price. 
The calculations shown in this paper were in part 
performed on the JUMP system of the Hochstleis- 
tungsrechenzentrum Julich. 

It is a pleasure to thank Joe Monaghan and Andrew 
Melatos for their hospitality during the last stages of 
writing this paper. 
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